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AB81A&CT 

TJs«  lBt«t;r&t«d  Ship  SbbIcii  Sy«t«oi  (tSDS).  is 
in  dnvttiopMCBt  to  pnnnit  pr«llninnry  sbip  dnnlse 
to  bo  accoiopliilind  using  Intsrnctlvo  computing. 
Curmntly  ISOS  donn  not  nddrons  nbip  dyncBicn  and 
tbis  rtfiort  nxploms  prograoM  in  tbin  dlscipUnn 
tor  inclusion  in  ISOS. 

A litorsturo  nnxrcb  «ss  psrtorasd  to  Idnntify, 
clnsslfy.  nod  tnbulntn  coig»ut«r  programs  involving 
sbip  dynamics.  The  search  included  programs  address- 
ing motion  prediction,  added  resistance  in  eaves, 
slamming,  deck  wetness,  and  maneuvering  for  dis- 
placQSiettt  hulls,  hydrofoils,  surface  effect  ships, 
axr  cushion  vehicles,  planing  bulls,  and  catamarans. 
Candidate  prograa&s  were  analyzed  with  respect  to 
their  input  data  requirements  and  the  relationship 
of  these  requirements  to  the  flow  of  data  among 
existing  design  programs  of  the  Integrated  Ship 
Design  System  (ISOS),  and  two  dynamics  programs 
for  displacement  hulls  are  reconnended  for 
integration  into  ISDS. 
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A«  NEED  FOR  THE  STUDY 

The  M&vy*s  Coaputer  Aided  Ship  DeslRitt  ftdd  Coimtructloii 
(CASDAC)  prograa  eddresites  the  t&eh  of  eaployiiw  coaputer 
progruus  to  iaprove  dealgn^  reduce  design  costs,  end  speed 
up  the  design  process.  CASOM;  is  generslly  concerned  with 
the  status  of  coaputer  progriaw  And  with  fostering  the  use 
cl  good  progreae  in  the  ship  design  process. 

A major  coaponent  of  CASXiAC  Is  the  integrated  Ship  De- 
sign System  (ISOS),  a coi'tecnion  of  coordinated  naval  design 
programs  which  interface  with  a coanon  data  base  to  form  a 
powerful  design  tool.  The  system  reduces  calendar  time  for 
the  design  process  by  facilitating  programs’  input  and  out- 
put and  by  use  of  Interactive  graphics  and  teletype  modes 
of  design  program  operation.  There  are  currently  no  ISDS 
programs  which  address  ship  dynamics.  There  Is  increasing 
interest  in  the  Navy  to  consider  the  maauvefing  and  sea°- 
keep  .ng  characteristics  of  ships  early  enough  in  the  design 
cycle  that  they  may  effectively  impact  the  design.  Inasmuch 
as  ISDS  does  not  presently  address  ship  dynamics,  it  is  pro- 
posed that  compatible  dynamics  programs  be  obtained  or  devel 
oped  for  inclusion  In  ISOS. 


B.  OBJECTIVE 

The  objective  of  this  task  was  to  determine  what  ship 
dynamics  programs  are  available  and  to  identify  any  caudi 
date  programs  suitable  for  inclusion  in  ISDS. 


C.  SUNNARY  OF  TASK 

The  task  perfomed  ess  to  sake  a literature  search  for 
cood>uter  program  Involving  ship  dynamics,  tahulate  the 
prograiis  found  froai  this  literature  search,  and  assess 
the  potential  usage  of  these  programs  in  ISOS. 


D.  SCOPE  OF  STUOY 

1.  Ship  Types  - The  literature  search  was  directed  toward 
finding  coagmter  programs  involving  the  following  vehicle 
types : 

Naval  and  Commercial  Displacement  Ships 
Hydrofoils 

Surface  Effects  Ships  (SES) 

Air  Cushion  Vehicles  .(ACV) 

Planing  Hulls 
Catamarans 

2.  Topics  of  Interest  -o  The  specific  topics  of  interest 
in  this  literature  search  are  the  following; 

notions  Predictions 
Added  Resistance  in  Waves 
Slamming 
Deck  Wetness 
Haneuvering 


. nTirrfilirh Hi'-ntiTTrnfrr 


11.  GENERAL  APPI^OACH 


A.  LITERATURE  SEARCH 

The  followlniE  literature  searches  were  coaducted  to  de- 
termine what  prograoM  are  available: 

A literature  search  of  the  DTNSBDC  library 

A Uarltlme  Research  Inforaation  Service  (MRIS) 

A Defense  Documentation  Center  (DDC)  literature 
search 

A Ship  Automatic  Retrieval  Project  (SHARP) 
literature  search 

In  order  to  assess  the  ship  dynamics  programs.  Informa- 
tion was  tabulated  on  each  program  found  in  the  above  litera- 
ture searches. 


B.  CONSULTATION  WITH  EXPERTS 


Experts  in  the  field  of  ship  dynamics  at  DTNSRDC  were 
consulted  to  determine  if  any  significant  ship  dynamics 
programs  were  missed  in  the  literature  search,  and  their 
recommendations  were  obtained  for  programs  to  be  integrated 
into  ISOS.  The  experts  consulted  were: 


Code 

Dr. 

Michel  Kazuu  Ochi 

1506 

Daniel  S.  Cieslowski 

1507 

Dr. 

Hanley  St.  Denis 

1513 

Dr. 

Wen  Chin  Lin 

1524 

Dr. 

Nils  Salvesen 

1552 

Geo  f f rey  G . Cox 

1568 

Nathan  K.  Bales 

1568 

Dr . 

Choung  Mook  Lee 

1572 

Sp«^c  lal  ist  - Naval  Architect 

High  Performance  Vehicles 
Program  Office 

Spec ial 1st  - Naval  Architect 

Design  Evaluation  Branch 

Hydrodynamics  Branch 

Surface  Ship  Dynamics  Branch 

Surface  Ship  Dynamics  Branch 

High  Performance  Craft 
Dynamics  Branch 
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III.  FINDIII6S 


A.  TABUUR  FINDIMGS  OF  LITERATURE  SEARCH 

The  fiadlnss  of  the  literature  search  are  presented  as 
two  tables  in  Appendix  B.  Table  1 is  a list  of  the  ship 
dynamics  cooiputer  programs  and  a brief  summary  of  each. 
Programs  listed  in  Table  1 are  those  for  which  the  library 
search  yielded  docusmntation  or  detailed  discussion;  Table 
I does  not  list  programs  for  which  the  search  produced  ref> 
erences  but  no  further  discussion.  The  entries  in  Table  1 
are  ordered  according  to  ship  type  and  topic  of  interest . 
with  the  newest  programs  in  each  grouping  listed  first. 

Table  11  is  a list  of  publications  encountered  in  the 
literature  search  which  reference  computer  programs  in  ship 
dynamics.  Including  some  of  the  programs  listed  in  Table  1. 
Table  11  is  also  ordered  according  to  ship  type  and  topic 
of  interest. 


B.  DISCUSSION  OF  FINDINGS 

1.  Six  Degree  of  Freedom  Programs  - For  conventional  naval  and 
commercial  displacement  ships,  there  are  three  6-degree  of 
freedom  (pitch,  heave,  roll,  sway,  yaw,  surge)  programs  which 
address  the  topic  of  motion  prediction: 

•The  OTNSROC  Ship  Motion  and  Sea  Load  Computer  Program^ 

•Program  SCORES  - Ship  Structural  Response  In  Vaves^ 

•The  MIT  6-Degree  of  Freedom  Program ^ 

All  of  these  programs  are  based  upon  the  compromise  theory.** 

•Superscripts  refer  to  references  in  Section  VII. 

•^Appendix  A provides  a brief  discussion  of  several  accepted 
theories  which  form  the  basis  for  various  dynamics  programs 


All  three  G-deffre^-of-freedon  progremiB  coaipute  the 
anplttude  utd  phase  for  the  pitch,  heave,  roll.  sway.  ynw. 
and  surge  notions,  the  vertical  and  horizontal  shear  forces, 
and  the  bending  and  torsional  ottaents.  Various  of  thev- 
prograns  also  Include  other  design  features  such  as  compu- 
tation of  hydrodynamic  pressure  distributions,  deck  wetness, 
added  resistance  and  slaamiog. 

The  6-degree  of  freedom  programs  require  as  input  a physical 
description  of  the  ship  including  hull  offsets,  ship  speed, 
ship  heading  relative  to  waves,  mass  of  the  ship,  mass  dis- 
tribution. and  a representation  of  the  sea  cordition. 

2.  Pitch  - Heave  Programs  ^ Ship  motion  programs  fur  the 
prediction  of  pitch  and  heav'  are  numerous*,  however,  most  of 
these  programs  are  based  on  one  of  two  theories  the  com- 
promise thfKiry  or  tho  Korvin-ICroukovsky  theory.^ 

The  Frank  Close-Fit  Ship  Hot  ion  Computer  Program^  com- 
putes pitch  and  heave  motions  for  ships  in  regular  and  irre- 
gular head  waves.  The  regular  wave  response;^  are  computed 
according  to  the  comprvMiuse  theory.  The  irregular  seaway  is 
represented  by  the  Pierson -huskowitz  Sea  Spectrum  and  statis- 
tical response  paraateters  are  computed  for  pitch  and  heave  as 
well  as  for  absolute  and  relative  vertical  displacement  veloc- 
ity and  acceleration  at  any  point  along  the  length  oi  the  ship. 

In  the  Frank  CIuse-Fit  Ship  Motion  Program,  there  art>  two 
methods  used  for  computing  the  two-dimensional  added  mass 
and  damping  coefficients,  the  Lewis-form  method  and  the  close 
fit  method.  The  Lewis-form  method  is  fast  but  it  is  only 
accurate  for  regular  (no  knuckles,  bulbous  bows,  etc.)  common 
ship  forms,  while  the  close-fit  method  is  time  consuming 
taut  accurate  for  any  station  shape.  The  user  has  the  optievn 
oi  selecting  the  best  method  for  each  station. 

•Francis  Ogilvie  and  Robert  Beck  present  a review  of  the  theory 
for  predicting  ship  motions'* 
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I.V  A liii  lil 


4ctothor  emuapie  of  tis«  pvourtta  itii 

“Cooputer  Aiidod  Prediction,  of  SeaEMpins  Perforscamee  la  Sl&ip 
Design"  by  T.  A.  Loukakin.  Tbe  theory  lised  la  thia  progna 
Is  essentially  tbe  strip  tbeory  of  Xorvia~Ki^iilu>vslcy  and 
Jacobs,  with  tbe  changes  presented  by  Gemtnea  and 
Beukeluan^.  The  program  will  calculate  the  response  oper- 
ators (the  ratio  of  the  ship  response  to  the  wave  asplitidt) 
of  pitching  and  heaving  motions  in  regular  waves,  it  will 
calculate  motions,  velocities,  acclerations.  relative  sotiona. 
relative  veir^cities.  bending  moments  and  shearing  forces  for 
any  position  along  tbe  ship,  and  it  will  calculate  the  cssan 
added  resistance  of  the  hull  form.  The  regular  wave  results 
are  used  in  conjunction  with  the  superposition  principle  to 
predict  iteakeeping  behavior  in  a seaway,  i.e..  in  a confused 
pattern  of  waves  of  various  sizes,  shapes,  and  directions. 

The  results  of  tbe  ship  behavior  in  a seaway  are  then  used 
In  conjunction  with  statistical  tbeory  to  predict  additional 
seakeeping  characteristics  of  a ship  such  as  probllity  of 
occurrence  of  deck  wetness,  propeller  racing,  and  slasmiag. 
Finally,  the  program  can  calculate  tbe  calm  water  resistance 
and  propulsive  coefficient  of  a hull  form  on  the  basts  of 
Series  60  or  Taylor  Series,  and  it  can  calculate  the  calm 
water  shearing  forces  and  bending  moaients.  In  general  the 
prediction  of  the  additional  seakeeping  characteristics— 
deck  wetness,  propeller  racing,  and  slamming — are  not  eon- 
sistance  with  experimental  results.  The  discrepancies  esn 
be  attributed  to  the  followAng  reasons: 
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a)  tiHi  dii^  at  a diiffaraat  dyaanic  draft 

uadar  ayaad  (aquat).  aad 
l»)  local  coadltioaa  pr«vallia«  at  tiM  boa  aad 
tba  atara  ara  aot  pradictad  by  theory . 

Aaotbar  aatbod  for  pradlctiac  pitch  and  haava  ia  ttia 
“ataadard  aarlaa”  approaeb.  Tba  tacbalqua  for  a staadard 
aarias  approach  ia  to  gaaarate  a faaily  of  bull  foraw  aith 
a ayataiaiitic  variatioa  of  all  seoaiatric  variablaa  ahich  ara 
axpactad  to  stroasly  iafluaaca  aaakaapia«.  Tha  parforaaaca 
of  ttaaaa  bulla  is  astabliabad  (by  nodal  testa,  full  scale 
taeaauraneat.  etc.)  &mA  data  are  recorded  to  produce  an 
eaienaive  data  base  ahich  nay  be  aubaequeotly  interpolated 
to  p«N»«lict  per fonaance  of  later  deataaa. 

The  Cctry  miA  Solution  l^rogran  uses  a standard 

aeries  approach  to  calculate  the  heave,  pitch,  relative 
notion  aad  absolute  acceleration.  As  ship  notions  do  not 
appear  to  be  aensittve  to  local  details  of  hull  shape,  it  is 
possible  to  select  s siiaplified  fanily  of  lifiathesuitical  foma. 
Each  natiunaatical  nodel  represents  all  hull  fonos  having  the 
aaiie  waterplaac.  profile,  and  sectional  area  curve.  The  aath- 
enatical  neodel  is  further  sinplified  by  adopting  a siaple 
syatenatic  set  of  eaterplanes  > each  baa  a parallel  nldship 
segwent  together  with  fourth  degree  ssolynonlals  which  swing 
into  the  centerline  at  the  fore  and  aft  perpendiculars 
without  a discontinuity  in  curvature.  The  profile  is 
rectangular,  so  all  sections  have  the  sane  depth. 
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3.  Eoll  Kotioc  - The  calcula.txon  of  the  force-aotlon  coef- 
ficicmttt  relating  to  roll  notion  Is  prinarlly  a nonlinear 
problem. 


Few  prograns  ealst  expressly  for  the  calculation  of  roll 
motion.  The  6 degree  of  freedom  programs  Include  roll  motions 
predictions  and  the  portion  of  the  program  calculating  the 
roll  motions  could  be  extracted  from  them.  The  6 degree  of 
freedom  program  by  Salveaen  et  al.  uses  an  inviscid  strip 
theory  to  coogiute  the  linear  components  of  wave  damping.  A 
quasi  linear  correction  term  was  then  added  to  account  for 
the  nonlinear  viscous  damping.  The  correction  term  depends 
on  frequency,  hull  geometry,  bilgekeel  dimensions  and  vis- 
cosity. A trj.al  and  error  procedure  is  used  to  solve  the 
equations. 
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Tb«re  Is  also  is  uss  at  OfUSBDC  tiie  Conally  Eoll  Prctstnui.^^ 
so  uopubllsheii  roll  In  iriilcb  lloosr  tlioory  is  uamnI 

to  reproMot  rollloc  Motion.  This  pro^rsM  gives  the  rolling 
motion  sod  moment  of  the  ship  In  regular  eaves  of  any  selected 
length  and  height,  and  this  may  be  combined  eitb  any  energy 
spectnim  of  save  slope  to  give  a spectrum  of  motion  la  irregu* 
lar  waves.  In  Irregular  eaves  the  response  amplitude  operators 
are  calculated,  and  the  ElfS  roll  and  mcment  are  gives. 

4.  Added  Eesistance  - A number  of  progrsms  exist  for  the  pre* 
dictlon  of  added  resist asce  In  eaves;  Appsndix  A describes 
the  theories  upon  which  four  of  these  progrsms  sre  founded. 

Since  added  renlstance  prograsMi  require  input  data  describing 
the  motion  of  the  ship,  the  programs  are  often  setup  to  run  in 
tandem  with  a ship  motion  program. 

The  added  resistance  programs  calculste  the  added  resis- 
tance in  regular  bead  seas  for  a number  of  wave  encounter 
frequencies.  Average  added  resistance  of  the  ship  proceeding 
at  a certain  speed  in  irregular  head  seas  can  then  be  predicted 
by  superimposing  the  programs'  added  resistance  in  regular  seas 
with  the  energy  spectrum  of  the  irregular  sea. 

5.  Slamming  and  Deck  Vetness  - M.K.  Ochi  and  L.E.  Matter^^ 
present  s method  for  prediction  of  slamming  charncterlstics 
and  hull  responses  of  a ship  at  an  early  design  stage,  given 
input  data  describing  the  Mtion  of  the  ship  and  certain  bull 
form  coefficients  which  are  significant  In  slam  prediction. 

The  predictions  include: 

e frequency  of  occurrence  of  slam  impacts. 

* elapsed  time  before  the  next  severe  impact  will 
occur. 

® magnitude  of  impact  pressure. 

^ the  largest  (extreme  value)  impact  pressure  most 
likely  to  occur. 

^ the  extreme  pressure  for  design  cons iderst ion. 

4 ship  speed  at  which  bottom  plate  daunage  is  most 

likely  to  occur. 
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• skip  It—  trcn  mlm  dmas** 

• littltiac  mp00d  tor  wbich  aUM  Impact  ia  tolarabla 
tor  tlia  crov. 

^ aitatial  dlstriUitloa  ot  ijipact  ptmumrm  aa<t  ramltiag 
torca. 

• »ato  hull  sirdar  raatmaaaa  such  sa  ahipplas  atrtUM 
sad  sccalsratioa  to  alaa  Ijaiwct. 

This  MSthod  has  sot  haaa  c<w|»lataly  prognuBaed,  hoaavar, 
portlosa  ot  tha  isathod  dasllag  with  traquaacy  ot  slaa  Isfuict 
hsva  baas  lataersiad  Isto  ship  aotloa  prograMS  sa  daalga 
tastiaras. 

Tha  pradlctloa  of  dack  vataaoa  has  also  baaa  Istagratad 
iuto  ship  notion  prograaai  as  a daslgn  taature.  Tha  ralatlva 
notion  batwaan  tha  ship  and  tha  aatar  surtnca  Is  obtalnad 
iron  tha  ship  notion  progrsn.  Standard  corraction  factors 
ara  applied  to  adjust  for  flara  and  tha  at fact  of  tha  ship's 
standing  wave  pattern. 


6.  hanauvering  - For  conventional  dlsplacanant  ships  there 

are  naneuvaring  progrsns  such  as  *‘A  Digital  Conputar  Technique 

Pr<»diction  of  Standard  Maneuvers  of  Surface  Ships'*  by  J.  Stron - 
13 

tejsen  and  "Directional  Stability  and  Control  of  ships  in 
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aestricted  Channels"  by  Uaruzo  Eda  . These  prograns  require 
as  input  hydrodynanic  force  and  monent  derivatives  which  nust 
be  obtained  froa  expensive  captive  nodal  tasting  techniques. 
Captive  BKMlel  tests  are  performed  by  naans  of  test  facilities 
such  as  the  rotating  am.  oscillators,  and  the  planar  antion 
nechansira.  The  hydrodynamic  derivatives  depend  largely  upon 
tha  ship  geometry  and  design  and  the  derivatives  differ  signify 
cantly  from  one  hull  fora  to  another.  The  hydrodynamic  deriva~ 
lives  are  by  definition  partial  derivatives  of  a force  or  moatent 
witb  respect  to  one  or  more  of  the  motion  variable.  To  obtain 

these  derlvati^’es  it  is  necessary  to  force  the  model  to  execute 

various  prescribed  motious  aud  to  ideasure  the  forces  and  moamnts 

as  functions  of  the  different  motion  paraaieters. 


There  ere  iqieclel  purpose  progrsnuB  for  various  classes 
of  ships  which  do  sot  require  the  hydrodynanic  derivatives, 
but  which  are  United  in  application.  These  prograsut  predict 
outaeuverlng  as  a function  of  various  paraaeters,  which 
define  ship  and  rudder  geometry,  and  which  define  operating 
variables. 

The  relationships  used  in  these  programs  were  derived 
by  use  of  regression  analysis  techniques  and  the  results  of 
the  progranui  are  only  applicable  to  ships  similar  to  those 
in  the  data  base. 

7.  Hydrofoil.  S£S.  ACV.  Planing  Hulls  - Programs  for  the 
prediction  of  ship  dynanics  of  hydrofoils,  surface  effects 
ships,  air  cushion  vehicles,  and  planing  hulls  are  time 
domain  programs.  Time  domain  programs  use  simulation  tech- 
niques to  compute  the  dynamics  of  the  vehicle  during  a suc- 
cession of  incremental  time  slices,  thus  giving  a time 
history  of  the  dynamics.  The  simulation  of  the  dynamics  is 
usually  carried  out  by  an  analog  program.  It  is  usually  not 
possible  to  write  such  programs  that  run  in  real  time.  The 
ratio  of  simulation  time  to  real  time  sometime  becomes  as 
large  as  100  to  1.  An  important  factor  in  the  motions  of 
these  ships  is  the  auto  pilot,  the  ship’s  real  time  control 
mechanism  which  interprets  Inputs  from  motion  sensors  and 
generates  directives  to  the  ship's  active  control  devices. 

The  auto  pilot  must,  of  course,  be  considered  in  the  program 
for  simulation  of  the  ship  dynamics.  There  is  no  general 
analog  program  for  ship  dynamics;  one  or  more  new  programs 
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Must  be  written  for  each  ship.  Host  of  these  ealstlag  pro- 
grans  are  not  doctSMnted.  Sons  of  the  specific  prograas  which 
are  docunented  are  listed  in  Table  1 of  Appendix  B. 

8.  Catanarans  - For  the  prediction  of  ship  action  and  loads 

15 

of  catanarans  there  are  prograas  such  as  MIT's  "CAT-5" 

This  prograa  is  capable  of  predicting  notions,  velocities, 
and  accelerations  of  catanarans  aa>ving  in  waves  at  any 
heading  angle.  The  wave  induced  forces  and  nooMuats  on  the 
catanaran  cross-structure  are  also  conputed.  These  responses 
can  be  obtained  in  regular  waves  as  well  as  long-crested  or 
short -crested  randan  waves.  For  the  latter  cases,  the  oue- 
paraneter  Pieraon-Moskowltz  Have  Spectrun  and  the  two-para- 
nater  Bretschneider  spectrun  are  used.  These  spectra  are 
defined  in  Section  l.d.  of  Appendix  A. 

The  basic  input  to  the  progran  includes  the  overall 
disMsnsions  of  the  cataanran.  the  fom  and  apacing  of  the 
deati-hulls,  the  longitudinal  and  transverse  weight  dis- 
tributions, the  general  heading  angle  and  the  description 
of  the  sea  condition  by  wave  spectra.  The  output  of  "CAT-S" 
consists  of  the  five  degrees  of  notion  (surge  neglected)  and 
the  wave  induced  forces  and  aoawtnts  on  the  cross-structure. 

There  is  also  a ship  notions  progran  for  catanarans 
being  developed  at  DTNSBDC  which  is  based  on  "Motion  and 
Resistance  of  a Low  Vaterplane  Area  Cataautran"  by  P.  C.  Pien 
and  C.  M.  Lee.^^  This  method  predicts  characteristics  of 
notion  and  hydrodynamic  loads  of  either  conventional  cata- 
nutrans,  or  small  laterplane  Area  Twin  Uull  Ships  (SVATH). 

The  equations  of  notion  for  catanarans  are  derived  in  the 
frequency  domain  under  an  assumption  of  a linear  excitation- 
response  relationship.  The  hydrodynamic  coefficients  Involved 
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in  tUe  eauntionn  of  Motions  arc  d«teniine<l  froa  strip  thsory. 
assuninR  slandor  gsoastry  of  «ach  bull  of  tb«  cataaaran.  Tbe 
effect  of  forvard  spcNtd  on  tbe  bydrodynaaics  coefficient  is 
treated  as  if  there  were  no  perturbation  on  tbe  fluid  due 
to  the  forward  aoveaent  of  the  ship.  An  apparent  underestiaa- 
tion  of  daaping  results  in  an  unrealistically  large  action 
aaplitude  at  tbe  resonant  frequency.  Tbe  introduction  of 
suppleaeatal  daaping  into  tbe  equation  of  aotion  is  needed  to 
achieve  a reasonable  prediction  of  eataaaran  notions.  The 
Suppleaeatal  daaping  is  introduced  into  the  equations  of 
aotion  so  as  to  be  linearly  proportional  to  tbe  oscillation 
velocity.  Prediction  of  statistical  averages  of  aotion 
aaplitude  for  the  cataafluran  in  irregular  seas  is  aade  by  using 
the  frequency  response  aaplitude  operator  in  conjunction  with 
the  Pierson-Moakowitz  sea  apectrua. 

The  probable  frequency  of  a water  contact  with  the  cross- 
deck structure  is  conputed  for  given  conditions  such  as 
significant  wave  height,  forward  speed,  and  height  of  tbe 
cross-deck  structure  above  the  waterline.  The  fomula  used 
is  based  on  a truncated  Rayleigh's  probability  distribution 
for  slaaaing  and  is  siailar  to  the  fomula  developed  by  Ochi 
(1964)  for  bow  slaanlog  of  aonohuil  ships. 

Respecting  added  resistance  for  cataaarans.  It  has  been 
shown  that  by  using  the  concept  of  an  "effective  bull  fora" 
the  added  resistance  for  cataaarans  can  be  predicted  with 
the  saae  aetbods  used  for  conventional  aonobulls. 
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IV.  SEOmEltATKM  FOR  OYH/litlCS  PROORAHS  IN  ISDS 


A.  RECOIINPIOATION  FROM  THE  OTNSROC  SHIF  PERFORMANCE  OEPARTfCNT 


Th«  Sblp  Perfonuuicv  0«partMittt  Jude  the  follcMPiug  r«c> 
onM0ttdatloiui  coac«riiiiifl;  programs  for  iatsfratioii  Into  ISOS: 

• Csadldste  programs  should  addrsss  pitch,  hsave  and 
loads  in  hsad  ssas  and  roll  in  bsam  s«as.  bscauss 
thsss  conditions  reprsssnt  the  worst  motions  and 
loads. 

* A pair  of  programs  > one  for  pitch-heave-loads  in 

head  seas  and  one  for  roll  in  beam  seas  - are  recosttend- 
ed  for  use  early  in  the  ISDS  design  process,  and  another 
similar  pair  for  more  exact  analysis  later  in  the  process 
These  four  programs  are: 

For  early  use: 

- For  pitch  and  heave,  the  Entry  and  Solution  Program 
which  uses  a standard  series  approach. 

- For  roll,  the  Conolly  Roll  Program. 


For  later,  wore  exact  analysis: 


- For  pi tcb-heave> loads,  the  Frank  Close  Fit  Ship~ 
Motion  Computer  Program. 


For  roll,  the  DTNSBDC  Sbip-Hotion  and  Sea-Load 
Computer  Program. 


These  recossmnded  programs  are  described  in  Appendix  C. 


B.  ISOS  DESIGN  PROCESS  CCNSIOERATIONS 

The  four  programs  recommended  by  the  Ship  Performance 
Departawnt  were  scrutlnzed  by  DTMSBDC  Code  1853  with  respect 
to  the  input  data  required  by  each.  The  sequence  of  ISDS 
programs  (See  Figure  1)  was  examined  to  Identify  the  earliest 
point  in  the  ISDS  process  at  which  data  is  available  to  allow 
execution  of  the  four  recommended  dynamics  programs. 


1 • »» lf/»a  wi 


Devlgn  Itachinery 
» 0tc . 


L. 


Dynamic  PrograouB 


1 
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Figure  1 > Portion  of  ISDS  Deaign  Process 
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Tils  fiadings  of  this  oaMulnatloa  are  as  follows; 


• The  principal  aignlflcast  difference  aaong  t2te  four 
prograsw  respecting  input  r<^4ulrefMnts  is  that  the 
hull  form  is  described  to  the  pair  of  earlf  programs 
merely  in  terms  of  basic  hull  dimensions  and  form 
coefficients,  whereas  the  two  mere  exact  programs 
require  a table  of  offamts  on  bull  stationa. 

• The  tiw)  more  exact  programs  also  roquire  input  Re- 
specting the  dynamic  properties  of  the  ship  (center 
of  gravity,  r^ii  .of  gyration,  etc.),  roll  damping 
and  bilge  keel  description.  Data  for  the  dynamic 
properties  and  roll  damping  can  be  satisfactorily 
estimated  from  known  properties  of  similar  ship 
types.  Feasible  bilge  keel  configurations  can  be 
assumed  and  verified.  All  these  data  can  be  com- 
piled by  the  engineer  subsequent  to  the  operation 
of  SCOl. 

• Hull  offsets  are  produced  by  the  second  ISOS  design 
procedure.  HULLCEN/BODYPLAN . These  offsets  are  not 
faired,  but  they  are  considered  to  be  acceptable  for 
use  in  the  two  more  exact  dynamics  programs. 


Given  that  the  two  more  exact  programs  can  be  exercised  so 
early  in  the  ISD8  design  process,  it  would  be  difficult  to 
Justify  the  effort  of  integrating  into  ISOS  all  four  dynamics 
programs.  However,  three  alternative  rationales  which  use  all 
four  programs  have  bees  identified  and  briefly  examined;  none 
of  these  rationales  are  considered  to  merit  further  development: 

(1)  Rationale:  The  two  early  prograaw  could  be  executed 

after  SCOl  and  in  parallel  with  or  prior  to  HULLGEN/BOPYPLAN . 
thus  verifying  the  seakeeping  characteristics  earlier  in  the 
design  process. 

Retort : UULLGEH/BODYPLAN  can  be  executed  in  1 - 3 

days  and  this  small  difference  does  not  .lustifv  the  extra 
implesHintation  cost. 

(2)  Rationale:  The  two  early  programs  could  be  incorporated 

within  SCOl  Itself,  and  seakeeping  considerations  would  be  avail- 
able early  enough  to  influence  the  selection  of  hull  dimensions 
and  form  coefficients. 
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Retort  Th«  two  early  procraBw.  aa  they  presently 
operate,  are  BULjor  procranw  in  their  own  rlRht.  and  require 
aodeat  CPU  tlM  to  execute.  Operation  of  these  prograns 
within  SCOl  to  exsalne  each  of  a large  aiatrlx  of  possible 
designs  wcwild  absorb  very  such  computer  time,  increasing 
aarhedly  the  operation  cost  of  SCOl.  Furthemore.  a aaajor 
coding  effort  would  be  required  to  join  the  prograss.  If 
Indeed  a measure  of  seakseplng  is  needed  la  this  formative 
design  step.  It  would  seem  that  a faster,  less  precise 
routine  should  be  developed  to  better  meet  the  need. 

(3)  Rationale:  At  the  stage  of  design  represented  by 

UULLGEN/BOOYPLAN  the  design  is  still  fluid,  and  the  dyna- 
mics programs  may  have  to  be  executed  many  times,  programs 
which  are  cheap  to  operate  are  desirable  for  the  early 
iterations. 

Retort : The  nu>re  exact  roll  program  runs  In  less 

than  one  minute  of  CPU  time,  so  there  Is  not  much  to  be  saved 
by  avoiding  its  use.  The  Frank  Close-Fit  Ship  Motion  Computer 
Program  may  be  operated  In  a "quick-and-dirty"  mode  utilizing 
a built-in  Lewis  form  to  represent  the  bull  form.  Execution 
time  would  presusmbly  be  reduced,  as  would  the  precision  of 
the  results,  but  the  user  would  enjoy  the  option  to  run  in 
either  mode. 

C.  RECOftlENDAnON  FOR  INTEGRAfED  SHIP  DESIGN  SYSTEM 

The  Frank  Close-Fit  Ship-Motion  Computer  Program  should 
be  used  to  investigate  pitch,  heave  and  loads  In  head  seas, 
and  the  NSRDC  Ship-Motion  and  Sea-Load  Computer  Program  should 
be  used  to  investigate  roll  In  beam  seas.  Both  programs  should 
be  Integrated  Into  ISDS  such  that  theyoould  be  executed  at  any 
time  after  the  ISDS  design  procedure  HULLGEN/BCOYPLAN  has  been 
run  (Figure  2). 
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V.  SiSfVm  AND  CONCLUSION 


A llteratiur#  ummxch  for  coi^putor  prosrus  lavolvlog  abip 
4yft«Miai  «ui  prforwd.  Tho  prograns  found  in  tbin  litornturo 
nnaxcli  am  liatad  with  a brief  description  of  each  in  Appendix  B. 

The  prograna  recoMsended  for  Integration  into  XSD8  are  as 
follows: 

*Tbe  Prank  Close-Fit  Ship  Motion  Coaputer  Prograw 
should  be  used  to  investigate  pitch,  heave  and 
loads  in  head  seas. 

*The  teSBDC  Ship-Motion  and  Sea-Load  Coa^^uter  Program 
should  be  used  to  investigate  roll  in  beam  seas. 

Both  prograaw  oiay  be  executed  very  early  in  the  ISDS 
design  sequence. 
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Presented  in  this  appendix  is  a brief  description  of 
sooie  widely  used  theories  which  have  been  progranmed  to  pre- 
dict ship  dynamics. 

1.  SHIP  MOTIONS 
a.  Strip  Theory 

In  strip  theory  the  ship  is  assumed  to  be  a slender  body, 
i.e. . the  beam  and  draft  are  small  cosMDared  with  the  length. 
Because  of  the  slenderness  of  the  body,  it  is  assumed  that 
the  major  fluid  motion  around  the  body  is  in  the  transverse 
directions,  there  is  not  enough  frontal  area  to  cause  impor- 
tant longitudinal  fluid  motion.  In  this  case,  the  flow  in  a 
transverse  cross-section  can  be  described  by  the  laws  of  fluid 
motion  in  two  dimensions. 

The  theories  which  are  presently  used  in  the  computer 
prediction  of  ship  motion  are  the  Compromise  Theory  and  the 
Korvin-Kroukovsky  Strip  Theory.  Both  of  these  are  strip 
theories  in  which  it  is  assumed  that  the  fluid  is  inviscid. 

The  viscous  effects  (damping)  are  small  and  can  be  disregarded 
except  in  the  case  of  roll  motion. 
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b.  CooproaiM  Tb«ory 

Tb«  comroalM  tlwory.  ditscribcd  by  SalvcuMi.  Ibek  and 

17 

fbitlttMa  r«sta  upon  tiu^  priacipftl  asstnptiou: 

• Tb«  bull  fora  Is  sssuflMd  to  be  s slender  body, 
i.e.,  the  beut  sod  draft  are  snail  coapared  with 
tbs  length.  More  precisely,  any  plaaes  tangent 

to  the  hull  should  be  nearly  paralled  to  the  lon- 
gitudinal sals. 

d It  Is  assuMd  that  there  are  no  significant  interactions 
beteeen  the  unsteady  floe  and  the  steady  disturbance 
caused  by  the  forvard  ant ion  of  the  ship. 

• It  Is  assuaed  that  the  velocity  potential  function  in 
the  final  equation  for  the  ship  notion  transfer  func- 
tions can  be  evaluated  atripwise. 

The  ccmpronise  theory  can  predict  the  heave,  pitch,  surge, 
ssay,  roll,  and  yae  notions  as  well  as  the  wave-induced  verti- 
cal and  horizontal  shear  forces,  bending  naawnts.  and  torsional 
nnnantfi  for  a ship  advancing  at  constant  speed  in  regular  waves. 


Under  ttui  assunptions  that  the  responses  are  linear  and 
haraonic  the  six  linear  coupled  differential  equations  of 
notion  can  be  written  in  the  following  fora: 

imrcla  forces  dasoiog  hydrostatic 
fores  fores 

where: 

Mjk  mre  the  coaponents  of  the  general ixed  aass  aatrix 
for  the  ship 

A.,  are  the  added  aass  coefficients 
jk 

are  the  daaplng  coefficients 
C..  are  the  hydrostatic  restoring  coefficients 

J ^ 

n^  are  the  translatory  and  angular  displacements 
h|j  are  translating  and  angular  velocities 

are  translating  and  angular  accelerations 

F.  are  the  complex  amplitudes  of  exciting  force  and  nostent 
«l 

u is  the  Ircqueucy  of  encounter 

i - vrr 

e ■ exponential  constant 


Tlie  six  couplsd  «qiuit;loiis  of  aotlon  roduco  to  two  lis- 
<tep«nckint  as  Mta  of  couplsd  iinoar  dlff  crest  lal  squat  loss 
with  frcquescy  asd  speed  depeadeat  coef flcieats.  This  re- 
duction Is  Bade  possible  by  recognixiac  the  lateral  syaaetry 
of  a ship.  One  set  of  three  coupled  equations  are  expressions 
for  surge,  heave,  and  pitch  uid  another  set  of  three  coupled 
equations  for  sway.  roll,  and  yaw.  Thus  for  a ship  with  lat- 
eral synaetry.  surge,  heave,  and  pitch  are  not  coupled  with 
sway.  roll,  and  yaw.  The  equations  for  the  wave  induced  loads 
are  expressed  in  terns  of  the  resulting  notions  and  the  de- 
rived hydrodynamic  coefficients. 

The  agreenent  between  the  conpromise  theory  and  experiment 
is  very  satisfactory  for  the  heave  and  pitch  notions  and  for 
the  vertical  loads,  in  oblique  and  lol lowing  waves  as  well  as 
in  head  waves.  Good  agreesMent  is  also  obtained  for  the  cou- 
pled sway  and  roll  notions  in  bean  waves.  Due  to  lack  of 
experinental  data,  it  has  not  been  possible  to  confirm  by 
experineut  the  tberoretical  notions  for  sway,  yaw,  and  roll 
in  oblique  waves.  Nevertheless,  the  good  agreenent  shown  for 
the.  horizontal  shear  forces,  bending  nonents  and  torsional 
atonents  in  oblique  waves  suggests  that  the  theory  may  also 
predict  tSie  horixoatal  notions  quite  well. 


c.  Korvia-Kroukovsky  Theory 

•hile  the  conpronise  theory  is  an  analytical  approach. 

g 

the  Korvin-Kroukovsky  theory  is  based  on  ingenious  and  care- 
ful physical  arguneate.  Because  the  Korvin-Kroukovsky  theory 
is  s strip  theory,  the  flow  around  the  ship  is  a 2-diiB8nsional 
problem.  Fi*om  the  solution  of  the  2-diaiensional  boundary  value 
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problMUi  of  tlMi  fIo«.  oB«  CAB  coMKmto  tlM  lijrdrodyBAaic  pros- 
Auro  dlAtrllMttlaA  oa  thm  afaip  caiummI  by  Ita  om  Botloa;  tlaiiA 
oaa  CAA  Alao  coapute  Addod  haaa  Aad  dimiiac  forcoe. 

It  la  oot  aaay  to  coaputa  tba  axcitiag  forcaa  daa  to 
WAvaa.  tmt  Xorvla-Kroukovaky  aada  a aiapla  Aaauaptloa:  la 

tba  AddadHtuMA  aad  daapiad  problaa,  ba  raplacaa  tba  valoclty 
And  AcoalaratloB  of  tba  ablp  croaa<-aaetloa  by  aa  affactlva 
ralAtiva  valoclty  aad  Accalaratioa  of  that  aactloa  altb  ra* 
apact  to  tba  InatAntaaaoua  vatar  aurfaca. 

FloAlly  it  vas  nacaaaary  to  iatrodaca  aomt  corractloa 
for  tba  affacta  of  forward  apaad.  Ooa  affact  la  tba  Dopplar 
ablft  la  tba  fraquaacy  of  eacouutar  of  the  ship  with  tba 
WAvaa.  Tba  aacoad  affact  la  that  tba  forward  not loo  of  tba 
ablp  cAuaaa  ao  apparaat  cross-flow  la  a traasvarsa  cross 
sactloa  if  tba  ablp  bappaas  at  saaw»  lostaat  to  hava  a aoa- 
zaro  pitch  aogla. 

Haaaoaably  accurate  predict ious  are  obtalaad  frosi  tba 
Korvio-Kroulcovsky  theory.  . The  theory  accurately  predicts 
the  heave  aad  pitch  not loos  as  well  as  the  vertical  shear 
forces  aad  beodlag  noautrats  for  regular  cruiser  stero  ships. 
Predict loas  are  less  accurate  for  ships  with  bulbous  bows, 
kouckles,  traoson  steros,  aad  other  irregular  features. 


d.  Motioos  la  Irregular  Seas 

A very  good  dascrlptioa  of  irregular  sea  conputatioas  is 
cootalaad  la  the  Fraafc  Cloaa-Flt  Sblp-Motloo  Coaputar  Progran 
by  W.  Fraak  aad  N.  Salvasao^  aad  is  quoted  as  follows: 

*'The  notloaa  of  tba  ablp  io  Irregular  seas  are  aaalyaad  la 

a atatlstlcal  naoaar  followlog  the  procedure  of  St.  Deals 
18 

aad  Plersoa.  Tba  statistical  approach  does  aot  give  the  actual 

tine  history  of  the  notioaa  but  rather  a statistical  dascrlptioa 
of  the  respoasa. 
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"Ttie  procedure  cun  be  outlined  in  short  ss  follows:  It 

is  sssuned  thst  both  the  excitation  (the  irregulsr  sen)  and 
the  i'esponses  (the  ship  ant  ions)  are  randon  processes  which 
can  be  assuned  to  be  stationary  (isussian.  Under  this  assusqF*- 
tion  an  adequate  statistical  knowledge  of  the  sea  elevation 
and  the  ship  responses  are  detemined  cosq>letely  by  their 
respective  enerfcy  epcctra.  The  energy  spectrrjB  is  a function 
which  specifies  the  fraction  of  the  total  energy  which  is 
associated  with  any  given  frequency  band.  The  sea  energy 
spectra  are  known  functions  furnished  by  the  oceanographer, 
while  the  response  spectra  must  be  computed. 

**The  second  fundamental  assumption  is  that  the  linear  super- 
position  principle  is  applicable  to  this  problem.  This  assump- 
tion lead.s  to  a simple  relationship  between  the  sea  spectra  and 

the  ship  response  spectrum.  More  specifically.  St.  Denis  and 
19 

Pierson  stated  that  'it  is  assumed  that  the  sum  of  the  res- 
ponses of  a ship  to  a number  of  simple  sine  waves  is  equal  to 
the  response  of  the  ship  to  the  sum  of  the  waves.*  Under  this 
assumption  the  ship  response  energy  spectrum  is 

3jj(w)  • (R<w)l^  . Sjj(w) 

where  R(u ) is  the  amplitude  of  the  response  to  a sinusoidal 
wave  of  unit  amplitude.  Sg(>^)  is  the  sea  energy  spectrum  and 
w is  the  wave  frequency.  Hence,  the  ship  response  spectrum 
S„(>^)  which  completely  describes  statistically  the  ship  motions 

n 

in  the  irreqular  sea  is  given  in  terms  of  the  response  ampli- 
tudes which  are  obtained  by  the  strip  theory  . . . and  the  sea 
spectrum  which  is  a known  fum  Ion.” 

The  following  three  sea  spectra  are  often  used  in  irregu- 
lar sea  calculations: 
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Mtfunaiui  tea  Sipactmn 
Th«  Bratseliattidar  tea  Spactnia 
the  Pittrsott^teakowltz  tea  Sjpectniii 


Tte  teuMUin  tea  SpectniM  la  a aiogla  paraaatar  apactrun 
which  daacrlbaa  fully  devalqpad  aaaa  aa  a functloa  of  wlad 
apaed: 

**  C« 

S^(w)  ••  aaa  a-^rcy  apactrun 

ahare  g ■ accalaratlon  of  gravity 

C • a conatant  (8.27  XIO*^) 

w " wind  apaad 


e « oxpotaatial  coaatant 

a • fraquaacy  of  «ava  pattam 

The  Bratachaeidar  ia  a two  paranatar  (aignlflcant  wave 
height  aad  oeau  period  of  waves)  aaa  apactrum  for  deacribing 
developing  and  decaying  aaaa 


S3(«) 


A“ 


4 


where  S^(a) 
A 


B 

hi 

T 

Tn 


u 

w 


sea  energy  spectrum 
4 v'h^VTta'* 

16»*/lto* 

significant  wave  height 

mean  period  of  waves 
frmiuency  of  wave  patteim 
hphed 


Ttitt  Pi«riK>ii»ilo«koirltz  S«a  Sp^ctniM  la  & aiacle  para«atar 
(alcalf leant  wave  baicht)  aea  apectniw  for  deacrlblag  fully 
developed  aeaa: 

where:  S,  (•)  - 4-  e ®* 


A 

B 

e.  Roll  Motion 


- .oosyi 


33.56/h:^ 

r 


The  linear  Methods  used  in  coe^>uter  prograas  for  the 
coaproalse  and  Korvin>Kruukovaky  theories  do  not  acceptably 
predict  roll  Motion  due  to  its  characteristic  non-linear 
nature. 


Although  much  work  has  been  done  in  other  fields  to 
develop  methods  for  the  efficient  solution  of  nonlinear 
problems,  little  use  has  been  made  of  these  methods  in  the 
field  of  ship  hydrodynamics. 

In  pure  roll  motion  there  are  two  principle  sources  of 
nonlinearity.  The  first  is  the  nonlinear  restoring  force. 

At  large  roll  angles  the  righting  moment  curve  may  not  remain 
linear  with  heel  angle.  The  degree  of  the  nonlinearity  de- 
pends on  the  ship  form.  The  other  source  of  nonlinearity  in 
roll  motion  is  the  damping  force.  Roll  damping  comes  from 
three  sources:  frictional  damping,  energy  lost  in  eddy  making, 
and  the  damping  due  to  wave  making.  Frictional  and  eddy  mak- 
ing damping  are  both  due  to  viscous  effects  and  are  essentially 
nonlinear  with  respect  to  roll  velocity.  Vave  damping  is 
associated  with  the  energy  carried  away  from  the  rolling  ship 
by  the  generated  waves  and  can  be  predicted  fairly  well  using 
linear  theory. 


In  the  technique  for  prediction  of  roll  used  in  the  Compro- 
mise Theory  as  described  by  Salveson,  Tuck  and  Fait  insen.  the 
viscous  rol l~damping  effects  are  Introduced  into  the  equation 
of  motion  as  a quasillnear  term.  This  term  is  the  product 
between  the  maximum  roll  velocity  and  a constant  which  depends 
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cm  til*  fr«qitMicy,  t9i«  vl«co*ity.  tb*  bll«*-ktt*l  dlnttasioas . 
and  tb*  bull  g«oii*try.  tb*  aaalnw  roll  velocity  is  *sti~ 
a*t*d  bttfor*  tb*  aotioii*  *r*  ca«|mt*d.  If  th*  diff*r*ac* 
b*t«e*B  tb*  •atlMBtod  and  tb*  cou|iut*d  maximum  roll  velocity 
is  too  lars*,  a a*«  value  for  tb*  siaalMMi  roll  velocity  aust 
be  estimated  aad  tb*  aotioa  tb*a  r*co«|>utttd  uatil  the  too 
COB verge. 


A liaear  theory  for  tb*  predict ioa  of  roll  is  presented 
by  J.  E.  Coaolly^^(1969).  la  this  theory  it  is  assuaed: 

« the  ship  is  wall-sided  at  the  waterline 

« the  buoyancy  forces  are  distributed  along  the 
length,  their  magnitude  depending  upon  the 
local  sectional  area  of  the  ship 

• the  local  force  is  assumed  to  act  normal  to  the 
local  wave  surface 

Integrating  over  the  length  of  the  ship  leads  to  an  effective 
wave  which  is  less  steep  than  the  true  wave,  and  which  has  a 
frequency  equal  to  that  with  which  the  true  wave  is  encountered. 

In  view  of  the  above  simplifying  assumptions,  it  is  consid- 
ered Justified  for  ease  of  solution  to  replace  the  ship  by  a 
siaqile  "equivalent**  hull  having  rectangular  sections  and  a 
constant  draft.  For  fine  warship  forms  the  design  waterline 
can  be  reasonably  represented  by  a parabola  having  Bttxiwuiti 
width  at  amldship  and  zero  half-breadths  at  the  ends.  Ellip- 
tical and  rectangular  waterplanes  forms  are  also  available. 

The  vertical  center  of  gravity  of  the  equivalent  bull  is 
placed  between  the  center  of  buoyancy  and  the  metacenter  in 
the  same  relative  position  as  that  of  the  real  ship  (l.e., 

18/^  is  the  same  for  both). 


The  "response  aaplitude  operator"  relates  the  roll 
responses  of  the  ship  to  various  regular  waves,  and  is 
defined  as  the  Maainuii  roll  angle  divided  by  wave  slope, 
have  slope  is  represented  as  a function  of  frequency  of 
wave  pattern,  wave  direction  relative  to  the  ship,  and 
ship  speed.  The  following  infomation  is  required  to  con* 
pute  the  roll  response  anplitude  operator  for  a particular 
ship  in  a particular  wave  pattern: 

•The  natural  period  of  roll  in  still  water 

•The  ratio  of  added  hydrodynamic  inertia  to 
the  ahlps  omumi 

•The  distances  between  vertical  center  of 
buoyancy  vertical  center  of  gravity  and 
laetaccnter 

•The  danping  monent  per  unit  roll  velocity 
as  a function  of  ship  speed.  This  moment 
can  be  expressed  non- dimensionally  in  terms 
of  a coefficient  of  roll  decay  which  defines 
the  manner  in  which  rolling  oscillation  would 
decay  after  removal  of  the  rolling  moment. 
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2.  AIQDSD  RESISTANCE 


Rempttctlng  tbe  topic  of  &dd«d  roslstiuice  in  vnvoa  for  con- 
vent lonnl  nhipe.  there  nre  four  distinct  theories: 

Potential  flow  solution,  by  Haruo 

Drift  force  approach,  by  Joosen 

Pressure  distributions,  by  Boese 

Radiated  euersy  approach,  by  Gerritawa  k Beukelaan 

These  theories  are  presented  and  coepared  by  Jorgen  Strom-Tejsen , 
20 

et.  al. 

The  theory  developed  by  Maruo  is  a potential  flow  solution 
to  the  problen^ in  which  the  velocity  potential  is  divided  into 
three  components  which  are  considered  separately.  A singular- 
ity distribution  is  used  to  represent  the  bull  form,  and  tbe 
wave  field  potential  consists  of  the  potential  associated  with 
tbe  regular  wave  field  and  velocity  potential  of  tbe  waves 
produced  by  tbe  singularities.  Haruo  has  shown  that  tbe  pitch- 
ing and  heaving  notion  of  tbe  ship  dominates  tbe  effects  of 
surge,  and  he  was  therefore  able  to  neglect  the  effects  of 
speed  variation. 

In  application  of  tbe  theory.  Haruo  has  further  approxi- 
mated tbe  ship  hull  by  a line  distribution  of  sources,  thereby 
allowing  the  evaluation  of  the  added  resistance  in  terms  of 
the  geometric  characteristics  of  tbe  ship  (described  at  each 
section)  and  tbe  measured  or  cooMiuted  ship  motion. 

The  added  resistance  theory  formulated  by  Joosen  is  based 
on  the  extension  of  a drifting-force-in-waves  theory  to  the 
case  of  a ship  under  way.  Tbe  drifting  force  problem  was 
completely  analyzed  by  Haruo  and  extended  by  Newman.  Haruo 
determined  tbe  drifting  force  as  a function  of  tbe  motion  of  tbe 
ship  and  the  wave  reflection  from  tbe  ship  hull.  Newman's 
contribution  was  an  extension  of  tbe  theory  to  Include  the 
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lott  oi  vertical  lao—at  and  tba  darivatioa  of  approxl- 
■ate  eapreasions  for  alandar  adlp  la  long  waves. 

Joosen  has  expanded  Maruo'a  drift  force  expression  in  an 
asymptotic  series  of  powers  of  a slenderness  parameter  and 
lias  truncated  the  series  after  the  first  order  terms.  The 
resultant  slender  body  approximation  for  the  drift  forces  is 
valid  for  short  waves.  The  theory  is  extended  to  account  for 
the  added  resistance  in  waves  at  ahead  speed  by  the  substitu- 
tion of  a wave  encounter  frequency  «diich  is  a function  of 
velocity  for  the  wave  frequency  appearing  in  the  analysis. 

Joosen  has  shown  that  the  drift  force  in  the  longitudinal 
direction  In  head  waves  is  dependent  only  upon  the  potential 
of  the  radiated  waves.  He  has  concluded  that  if  the  ship 
motion  amplitude  is  of  a similar  order  of  ouq^nitude  over  the 
entire  range  of  encounter  frequencies,  the  wave  diffraction 
effect  may  be  neglected  everywhere  except  for  very  ssmll  vimm 
where  it  beccwms  dominant. 

Joosen' s final  expressiona  are  equivalent  to  those  derived 
by  Havelock  except  for  the  inclusion  of  a wave-motioEi  coupling 
term.  The  added  resistance  is  determined  through  a knowledge 
of  the  ship's  frequency  dependent  dampinig  f icient»  as  well 
as  information  on  the  heaving  and  pitchinj^  motion  and 
relationships. 

Boese  has  recently  studied  the  added  resistance  problems 
using  an  approach  similar  to  the  original  work  perfonued  by 
Havelock,  but  considering  a more  sophisticated  approach  to  the 
Initial  motion  data.  Boese  states  that  since  the  three-dimen- 
sional pressure  distribution  on  a ship  in  waves  cannot  pre- 
sently be  determined,  it  is  necessary  to  employ  linear  strip 
theory  for  computing  the  ship  motions,  und  hence  also  for  the 
pressure  distribution. 


Tti«  pressure  Xorces  scting  oa  s ship  ere  divide4  by 
Boese  into  two  segaecta: 

• The  fluctuating  force  caused  by  the  wave  fif.id 
and  the  heaving  and  pitching  Motion  over  a fixed 
surface,  i.e. . the  hull  below  the  static  water- 
line. 

« A correction  to  the  fixed  surface  force  to  account 
for  forces  over  regions  of  the  hull  which  are  not 
continuously  iHBOiersed.  For  this  correction,  dynan- 
ic  forces  are  neglected  and  a linear  pressure  dis- 
tribution is  assumed  to  be  effective. 

Computing  the  tine  average  contribution  to  the  longitu- 
dinal force  costpletes  the  calculation  of  added  drag. 

Bcese  Justifies  his  method  on  the  ground  that  first  order 
effects  dominate  in  notion  responses. 

The  method  used  by  Gerritsma  and  Beuhelmau  to  calculate 
added  resistance  is  to  compute  the  energy  flux  radiated  from 
the  hull.  The  energy  is  found  in  terns  of  the  sectional 
damping  coefficient  and  the  vertical  relative  water  velocity 
between  the  ship  and  the  wave.  The  radiated  energy  is  then 
equated  to  the  work  being  done  by  the  ship.  Gerritsma  theory 
has  the  advantage  that  the  added  resistance  is  easy  to  compute 
if  one  is  already  coogfuting  the  notions  fron  strip  theory. 

Also  the  radiated  energy  approach  gives  perhaps  the  easiest 
description  of  resistance  in  quartering  and  bean  seas.  The 
added  resistance  in  oblique  seas  is  found  by  considering  the 
energy  radiated  fron  each  side  of  the  ship  separately. 
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3.  SLAMMING 


Mlcbel  K.  Ochl  &ad  Levis  E.  Hotter  present  s nw»tbod  to 
predict  the  slsaniOE  characteristics  and  hull  responses  of 
a ship  at  an  early  design  stage.  A statistical  approach  to 
the  developoMnt  of  a prediction  ate t hod  is  taken. 

The  following  is  a brief  outline  of  the  overall  predic- 
tion aetbod.  Two  prerequisites  have  to  be  prepared  in  advance. 
One  is  the  description  of  ship  Motions  and  the  other  is  the 
evaluation  of  coefficient  associated  with  slaa  impacts.  A 
coaputer  program  to  evaluate  the  coefficient  for  a given 
section  shape  is  available  in  reference  30. 

With  these  two  prerequisites  in  hand,  the  following  pre- 
dictions can  be  made: 

a.  Frequency  of  Occurrence  of  Slam  Impact.  This  can  be 
easily  obtained  from  the  inforaation  on  ship  motion  and 
velocity  relative  to  waves. 

b.  Limiting  Speed  for  Which  Slaaaing  lapact  is  Tolerable. 
From  results  obtained  in  Item  a as  a function  of  ship  speed, 
the  limiting  speed  for  which  slaa  iapact  is  tolerable  - here- 
after referred  to  as  tolerable  speed  for  slaa  impact  - is 
estimated  as  the  speed  for  which  either  the  probability  of  slaa 
inqiact  at  Station  3 reaches  a level  of  0.03  or  the  significant 
amplitude  of  the  vertical  bow  acceleration  reaches  a level  of 
0.4  g.  This  subject  is  discussed  later  in  the  section  on  ship 
speed  and  slamming  damage. 

c.  Time  Interval  Between  Successive  lapacts.  From  results 
obtained  in  Item  a the  time  interval  between  successive  impacts 
can  be  estimated. 
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d.  Sl—aittg  Prammr*.  thm  firobftbility  function  necesn&ry 
for  predict lac  lapact  presaure  la  eatabllabad  from  the  ahip 
Botloea  aad  coefficient  aaaoclatad  with  alan  ispacta  thua  the 
average  aad  algnlf leant  valuea  of  Impact  preaaure  can  be  pre- 
dicted. 


e.  Elapaed  Time  Before  the  Heat  Severe  lapact.  The  time 
between  one  severe  Impact  and  the  next  at  a specific  location 
along  the  ship  length  can  be  estimated  from  the  results  obtain- 
ed In  Items  a and  d. 


f.  Most  Probable  Extreme  Pressure.  By  applying  order 
statistics,  the  magnitude  of  the  largest  impact  pressure  most 
likely  to  occur  in  a specified  ahip  operation  time  in  a given 
sea  can  be  predicted  at  any  location  along  the  ship  length  from 
the  results  obtained  in  Item  d. 


Ship  Speed  at  Which  Bottom  Plate  Damage  is  Most  Likel: 


to  Occur.  By  equating  the  pressure  laagnitude  obtained  in  Item 
f to  that  which  will  cause  permanent  set  of  a rectangular  plate, 
the  speed  is  estimated  at  which  bottom  plate  damage  is  most 
likely  to  occur. 


h.  Extreme  Pressure  for  Design  Consideration.  The  magni- 
tude of  extreme  pressure  for  design  consideration  can  be  esti- 
mated by  applying  order  statistics  to  the  probability  function 
obtained  in  Item  d.  The  extreme  value  is  controlled  by  a pre- 
assigned small  probability  of  being  exceeded,  which  the  designer 
specifies. 


i.  Ship  Speed  Free  From  Slam  Damage.  By  equating  the 
pressure  magnitude  obtained  in  Item  h to  that  which  will  cause 
permanent  set  of  a rectangular  plate,  the  attainable  (maximum) 
speed  below  which  no  plate  damage  would  occur  is  estimated. 

The  probability  of  occurrence  of  damage  is  the  small  number 
assigned  in  item  h. 


J.  Sl«»  liipact  Force.  The  aeffnltikile  of  Inpact  force 
cen  be  estlnated  fron  Itene  f or  h,  takiag  spat lei  distri- 
bution and  travellnc  tine  of  the  pressure  into  consideration. 
The  force  evaluated  using  the  extreaw  pressure  In  Item  h is 
used  for  design  consideration. 

k.  Main  Hull  Girder  Response  to  Itapact.  Using  the  in- 
fomation  obtained  In  Iten  <1  as  an  Input,  the  nain  hull 
girder  responses  such  as  whipping  stress  and  deceleration 
due  to  lopact  are  estlnated  by  solving  a aathenatical  node! 
representing  the  hull  structural  chara'^terlstlcs.  Available 
conputer  prograns  for  dynanlc  response  of  the  hull  may  be 
used  with  necessary  alterations. 

According  to  M.  K.  Ochl  and  L.  £.  Hotter  (1973).  a com- 
plete method  to  estimate  In  the  design  stage  the  slanning 
characteristics  and  hull  responses  of  a ship,  including  the 
estimation  of  ship  speed  for  which  damage  on  bottom  plating 
could  possibly  occur,  is  not  available  at  present.  There 
are.  however,  numerous  papers  on  methods  to  obtain  limited 
information  on  slamming  Including  both  theoretical  and 
experimental  approaches. 


4.  MANEUVERING 

In  order  tr  calculate  the  path  of  a ship  the  hydrodynamic 
force  and  moment  derivatives^ which  are  by  definition  partial 
derivatives  of  a force  or  moment  with  respect  to  one  or  more 
of  the  motion  parameters,  must  be  known.  The  only  way  to 
obtain  these  derivatives  is  from  captive  model  tests.  Captive 
model  tests  are  performed  by  means  of  test  facilities  such  as 
the  rotating  arm,  oscillators,  and  the  planar  motion  mechanism. 
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Tlie  hydrodyaanic  <torivatlv«aidapead  largaly  upoa  the  ship 
gttoiMtry  sad  dssica  sad  la  csaeral  they  differ  algal f least ly 
frcM  oae  hull  fora  to  aaother.  To  obtain  the  hydrodyasalc 
derivatives  it  la  aeceaaary  to  force  the  aiodei  to  execute 
various  prescribed  aotloas  and  to  Measure  the  forces  sad 
Monsata  as  fuactioas  of  the  dlffereat  aotioa  parsioeters. 

The  hydrodyasnlc  force  sad  MosMsat  derivatives  are  used  in 
coabinatlon  with  the  equations  of  notion  to  calculate  the  path 
of  a ship. 
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mmixi 

TABULAR  FlKDllUiS  OF  LITCRATUlt^  SEARCH 

i^p«fiellx  B la  praaaatad  aa  t«o  tablaa.  Table  X la  a list 
o£  tba  cca«»utar  protnuaa  deallag  wltb  ahip  dyaamlca  fouad  la 
tbe  lltaratiure  aaarch  aad  iacludaa  a brief  aiaiaiary  of  each 
progran.  Table  XI  la  a Hat  of  pubXlcatioaa  eacoiiatered  la 
the  literature  aaarch  which  refereace  coaputer  prograua  ia 
ship  dyaamioa. 
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COMPUTER  PROGRAMS  IN  SHIP  DYNAMICS 
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tamnY  AMP  sourriow  p»xauts  (esp) 

SSP  utiUxM  A siAttdard  ««rt«s  approach  to  eatlaiate  pitch 
aad  lioav*  awtloaa  in  head  seaa.  It  was  develop«^d  by  N.  K. 

■alM  of  OTMllDC. 

EdS^  alluwa  deteraiaatlon  of  trenda  In  aeakaepiac  responaes 
with  chaasea  in  hull  t«<iHetry  at  an  early  atac®  la  the  design 
proceaa.  The  neaheeplag  reaponaen  conaldered  are  heave,  pitch, 
and  the  aaplitude  and  acceleration  of  vertical  Motions  near 
the  bow.  The  user  inputs  the  aignlf leant  wave  height  (average 
of  the  one-third  highest  waves)  representative  of  the  sea 
state  of  interest,  and  the  prograa  uses  the  Plerson-Moskowi ts 
wave  spectnia  to  coaipute  each  response  at  the  Input  speed  and 
at  speeds  corresponding  to  Froude  nunbers  of  0.1.  0.2.  and 
0.3. 

ESP  requires  the  following  input; 

Ship  geoaietry: 

• length  between  perpendiculars 

• beaw 

• draft 

• nuwinal  sectional  area  coefficient  at  forward 

perpendicular 

• noainal  sectional  area  coefficient  at  aft  perpendicular 

• waterplane  coefficient  forward  of  midship 

• waterplane  coefficient  aft  of  midship 

« midship  section  coefficient 

Significant  wave  height 

Ship  speed  (knots) 


Set  of  response  weighting  factors  (optional); 

Vl-weighiiog  factor  for  ms  heave 

V2-welghting  factor  for  rms  pitch 

V3-weightlng  factor  for  rms  relative  motion 

Vd-weighting  factor  for  rms  absolute  acceleration 


0 


' fiiiwiwii-n  '| 
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TIm  output  froa  ESP  Is: 

SMS  heave  asplltude  (tt) 

SMS  pitch  aaplltude  (ft) 

SMS  relative  aotloa  (ft)  | at  10%  of  shl»  lessth 
SMS  acceleration  (ft/sec^))  (Station  2) 

Seakeeping  efficiency,  a coopoalte  oieasure  obtained  by  a4>plylng 
weighting  factors  to  the  above  four  notions. 

The  weighting  factors  are  used  only  to  eaphaslze  the  rela- 
tive Inportance  of  the  different  notion  responses  In  the  com- 
putation of  seakeeping  efficiency,  an  enperical  measure  of 
overall  seakeeping  ability.  The  assignment  of  weighting 
factors  Is  a value  Judgement  which  the  designer  must  make 
on  the  basis  of  the  mission  requirements  which  he  is  attempt- 
ing to  satisfy.  If,  for  instance  all  four  responses  are  con- 
sidered to  be  equally  istportant,  he  should  assign  a weight  of 
1 to  each  response.  If  on  the  other  hand,  relative  motion  is 
considered  to  be  significantly  more  important  than  the  other 
responses  he  might  assign  it  a weight  of  5 and  the  other 
responses  weights  of  1. 

The  current  data  base  is  limited  to  single  screw  cargo 
hulls  with  either  conventional  sterns  or  transom  sterns.  To 
extend  this  data  base,  say  to  destroyers,  would  cost  approxi- 
mately 60X  per  class  of  ships. 

On  the  CDC  6700  computer,  this  program  requires  71.200 
octal  words  to  load  and  a maximum  of  40  cpu  seconds  to  run, 
however  this  program  has  not  been  optimized  for  either  min- 
imum memory  or  running  time. 
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COWOLLY  aOLJL  tmoCttAM 

*thm  CoBolly  Boll  Prograoi  is  sn  uapubUsli«d  roll  procran  la 
wlUch  It  is  sasuasd  that  tha  ship  is  wall  siM  at  the  water- 
lias  and  that  distributed  aloac  its  leacth  are  buoyaacy  forces 
whose  sacoltudes  are  deterslaed  by  the  local  sectioaal  areas. 
The  local  buoyaacy  forces  are  assmatd  to  act  uorsal  to  the 
local  wave  surface.  This  theory  lavolves  a auaber  of  sl^ll- 
fylac  aasusptioas,  aad  la  view  of  these  it  is  coasidered 
Justified  for  ease  of  solutloa  to  replace  the  ship  by  a 
slaiple  equlvaleat  shape  havlas  rectasdular  sect  loss  a 
coastaat  draft.  For  flaw  warship  foraa.  the  waterllae  width 
caa  be  reasoaably  repreaeated  by  a parabola  haviag  — — 
half-breadth  at  the  mid-leacth  aad  aero  half-breadth  at  the 
eads.  This  shape  teads  to  over-estlsate  the  width  of  the 
forebody,  aad  uader-estlauite  that  of  the  afterbody.  Ellipti- 
cal aad  rectaasular  waterplaae  shapes  are  also  available. 

This  progras  gives  the  roll  lag  aotloa  aad  rlghtiag  am 
soseat  of  the  ship  la  regular  waves  of  aay  selected  leagth 
aad  height.  Aay  eaergy  spectrus  of  wave  slope  auiy  also  be 
specified  to  give  a distribution  of  uiotloo  la  irregular 
waves,  la  Irregular  waves  the  respoase  aiaplltude  operator, 
the  RMS  roll,  aad  the  heel  lag  soweat  are  calculated. 


The  laput  for  the  Coaolly  Kotl  Progras  is  as  follows: 


Ship  geoaurtry; 

•leagth  between  perpendiculars 

•beas 

•draft 

•center  of  buoyancy 
•setacentric  height 


•natural  roll  period 
•coefficient  of  roll  decay 

•ratio  of  hydrodynasic  inertia  to  ship  inertia 
Speed 


Significant  wave  height 
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Ttui  output  for  ttilu  prosraa  Is  tbs  BMS  Boll  Baplltuds 
of  tbs  ship  to  rsculsr  wsvss  for  ov«ry  IS  of  hssdisc- 

On  ths  COC  6700  tbs  progrsa  rsquirss  16.200  octnl  sorOs 
to  loud.  A typical  run  to  msnlyass  thrss  ship  spssds  st 
bsntfiJMSof  svsry  IS  dscrsss  rs^sulrM  about  25  ssconda  of 
coaputsr  tlas. 
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Mm  wmm  mtmtmM  la  r^racrapli  IV.  B. 2.  tiMrr  ara  a largi^  aua> 
bar  ot  pmcranai  tor  tba  caloulatloa  of  pitch,  boava  aad  loada. 
tiMMaa  rrinifTiMia  ara  baaad  oa  tba  aaaa  tlMori>-'*>.  and  all  calcalata 
approaiaataiy  tba  aaaa  outputs.  Tba  proitrao  iaput  tovmt,  alas, 
aad  ruaalac  tlaw  mm3>  vary  from  prograai  to  procrao  but  tba  calcu- 
latloaa  art  all  aladlar.  tba  apacltlc  probrao  recrwaiMiad  tor 
parfonilaff  tba  calculatloa  la  tba  7raak  Cloaa  Fit  Sblp  Motioa 
Coaiputar  Prosnua.  Tbia  procraia  la  racaamadad.  avaa  tbouBb  It 
la  difficult  to  say  tbat  oaa  proarao  ta  battar  tbaa  tba  otbar. 
bacauaa  of  fanlllarity  of  aavy  paraoaaal  «Ub  ibis  prograa. 

Tba  Fraab  Cloaa  Fit  Ship  Hotioa  Cosputar  Procrao  does  sot  at 
prasaat  calcuiata  aoaui  daaiffa  iafonaatioa  (aliwsini,  dock  aat> 
aaaa.  propallar  aoNirvaaca)  tbat  aould  ba  uaaful.  lacludiag 
thaaa  daalga  icaturaa  would  be  only  a aauin  extaasion  of  tba 
praaaat  progran. 

Tba  gaoaetry  of  the  aulp  ataticma  oay  ba  daflaad  by  two 
diffarast  natboda,  aaacly: 

a.  by  giviag  about  aigbt  offaat  polats  for  aacb  atatloa.or 

b.  by  giving  tba  baa».  tba  draft  aad  tba  araa  coafficlaat 
of  tba  cross  aactlons 

Note  tbat  it  la  eot  aacvsaary  that  tba  same  aatbod  ba  uaad 
for  the  autira  ship. 

If  a station  ia  rapraaantad  by  offaat  pointa,  tba  uaar  than 
baa  tba  fraadoui  to  aalect  betaean  the  cloaa<>flt  oatbod  >r  tba 
Lawia-fora  oetbod  in  tba  coagHitatioa  of  tba  aactioaal  addad 
aaaa  and  daapiag  cocf  ficiants.  Vban  tba  Lawia>fora  aatbod 
is  salectad.  tba  coaputar  scans  the  offsats  for  aacb  aactioa 
to  dateraiaa  the  local  bamt.  draft,  aad  aactioaal  araa.  and  a 
Lawis  fora  awtebiag  tboaa  paraaieters  is  salactad  for  tha  sac> 
tioo. 
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la  additioa  to  tha  lafonMttloa  about  tte  gaoHatry  ot 
tba  abip.  tba  uaar  alao  auat  aalaet  tba  coa41tlooa  for  ablcb 
ba  aaata  tba  coaputatloaa  parfonM4.  For  ragular  aava  out- 
put ba  auat  aalaet  tba  abip  <^paaga  aad  tba  raaga  of  aacouatar 
firaquauciaa. 

for  Irfagular  aaa  coosiutatioaa  it  ia  aacaaaary  to  apaci- 
fy  tba  abip  apaate.  abip  laagtb.  tba  aigalficaat  aava  baicbt, 
aad  tba  loagitudiaal  locatioaa  at  abicb  aotioa  raapoaeas  ure 
to  bo  calculatad. 

Yba  Piaraoa-lioaboaita  aaa  apactruu  la  aapraaaad  la  tbiu 
prograa  aa  a fuactioa  of  tba  aigaif leant  wave  balabt. 

Tba  prograu  la  daalgaad  ao  that  Irregular  aaa  reaponaea 
can  be  couputad  at  ala  dlffaraot  abip  leagtba  for  the  saue 
given  abip  foru  without  any  cbaagea  to  thoae  data  cards 
wbleh  give  tba  gaoautry  of  tba  ship  foru. 

The  laput  for  the  Fraab  Close  Fit  Ship  Motion  Coeputer 
Frograa  la  aa  follows: 

Ship  geeawtry: 

* loagitudiaal  radius  of  gyration 

* length  between  perpendiculars 

* geoaatrlc  deacrlptlon  of  an  adequate  aimber  (leas 
than  SO)  of  ship  stations,  either  as  offsets  or  as 
the  beau,  draft,  and  area  coefficient  for  each 
station 

For  regular  waves; 

* Froude  nuaiber 

* range  of  eacounter  frequencies 

For  irregular  waves: 

* significant  wave  heights 

* range  of  ship  speeds 

* points  along  ship  for  which  wot  ion  responses  are  to  be 
calculated 
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Output  froui  this  progruM: 


* coefflclunta  aad  forcas  In  tba  aquatloua  of  aotloa  for 
pltcb  attd  tiaava 

For  ragular  tMbvaa: 

* pitch  aad  haava  aaplltudaa  and  phaaaa  over  ranga  of 
ancountar  fraquancias 

* plota  of  pitch  and  haava  aaQ>litudaa  and  thalr  phaaa 
anglaa  aa  a function  of  wava  langth 


For  Irrajgular  wavaa:  ^ 

* tba  avaraca.  algnlf leant,  and  1/10  bighaat  aaplltudaa  , 

of  pitch  and  haava  | 

• tha  abaouta  and  ralativa  vertical  dlaplacaoumt  valocity  . 

and  accalaratlon  at  any  point  along  tha  langtb  of  tba  ] 

ahlp  g 


Propart laa  of  input  hull  form: 

• block  coafflclent 

• diaplacamant 

• longitudinal  cantar  of  buoyancy 

• longitudinal  contar  of  flotation 

• aactional  beam 

• sectional  draft 

• aectlraal  araa  coefficient 

• plot  of  section  shapes 

The  typical  running  tioM^  for  this  program  using  20  sta- 
tions and  40  frequencies  is  500  CPU  seconds  on  tha  CDC  6700 
besavar,  if  all  Lewis  form  sections  are  used  the  running  time 
Is  reduces  to  about  50  CPU  seconds.  The  running  time  is 
proportional  to  tha  number  of  stations  and  frequencies  used. 
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PTMSRDC  Shlp«Motlon  and  8«a  Load  Progrm 

This  prosrui  Is  s gsnsrsl  six  dsfrss  of  frssdcsi  aotios 
prognuB  froa  wblcli  tbs  portion  of  tbs  prograa  for  calculst> 
Ing  roll  csn  bs  rssovsd.  Tbs  rssson  for  not  using  this 
progrsst  as  is  for  tbs  calculation  of  pitch  and  bsavs  as 
ssll  as  roll  is  that  tbs  prograa  is  a rsssarch  tool  pro~ 
viding  mich  aors  Infonnation  than  is  nesdsd  for  dssign,  and  it 
is  vsry  tiae  consuaing  to  run. 

The  NSfiOC  Ship  Motion  and  Sea  Load  Coaputer  Prograa  in 
its  present  fora  predicts  the  roll  aotions  for  conventional 
cruiser  stem  bull  formm  at  moderate  speeds  in  beaa  seas 
with  reasonable  accuracy^. 

In  order  to  calculate  roll  aotions  which  are  reasonably 
accurate  in  the  frequmcy  range  close  to  resonance  it  is 
necessary  to  obtain  good  estlaates  for  the  viscous  roll 
damping  coefficients  used  in  the  equation  of  aotion.  The 
total  and  sectional  viscous  roll  damping  coefficients  are 
not  generally  available  for  most  ship  fonts  and  they  are  ex- 
tremely difficult  to  obtain. 

Non-linear  viscous  roll  damping  is  introduced  in  the 
six-degree  of  freedom  equations  of  motion  in  the  fora  of  a 
quasi- linear  term  which  is  added  to  the  linear  damping  co- 
efficient. This  quasi-linear  term  depends  on  tbs  frequency 
of  encounter,  the  viscosity,  the  bilge  keel  diaensions,  the 
hull  geometry,  and  the  aaxiimim  roll  angle. 

The  program  usee  a trial  and  error  procedure  for  solving 
roll  aotions.  First  the  roll  responses  are  computed  using 
an  estimated  value  for  the  maximum  roll  angle  keeping  the 
wave  slope  constant,  then  the  program  compares  the  maximum 
cosqjuted  roll  amplitude  with  the  estimated  value.  If  the 


dlff«r«iice  Is  Isrgsr  thsa  ose  <Ssgrs«,  s ss«  vslus  tor  tbo 
ssslMWMS  roll  angle  is  estlaated  by  the  prograa  sad  the  cob- 
putations  are  repeated  until  the  difference  between  eetiaated 
and  coaputed  BanlBUB  roll  anplitude  is  leas  than  one  degree. 

The  input  for  the  roll  portion  of  the  NSBDC  Ship  Motion 
and  Sea  Load  Conputer  Prograa  is  as  follows: 

• offsets  on  21  stations 

• mass  of  ship 

« roll  and  yaw  radii  of  gyration 

• aass  product  of  Inertia  about  the  X and  Z axes 
« vertical  and  longitudinal  centers  of  gravity 

• total  length  of  submerged  portion  of  the  hull 

• classification  of  each  section  shape  (e.g..  deep  V, 
extrcMuely  rounded  sections,  etc.) 

« dead  rise 

• bilge  radius 

• bilge  keel  description  (See  Figure  3 ) 

• length  and  width  of  bilge  keel 

• angles  a and  0 

• girth  distance  (G)  from  bilge  keel  to  DVL 

• distance  (D)  frooi  bilge  keel  to  a longitudinal  axis 
through  the  CG. 

The  output  from  this  program  is: 

• the  coefficients  and  forces  in  the  roll  equation  of 
motion 

« the  roll  amplitude  as  a function  of  frequency  of 
encounter 
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Figure  3 
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- Typical  Shi] 
GeoMtric  Pi 
the  Locatloi 
Keels 


• tbm  roll  TmpoBM  opermtora  (glv«s  roll 

Mplituiia  rnn  £ ftasctioii  of  wwvm  halgitt) 


• tha  roll  nooHlljMiuiioiuil  traasfor  functions  (givos 
roll  aaiplituiitt  &•»  & funetioa  of  w«v«  alope) 


Tfcne  •stlaatatl  tiaa  to  run  t&iia  pmgrsui  for  roll  notioo 
la  less  thss  60  sscosds. 
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INITIAL  DimiBIITIOM 


NAV8EA 

1 SEA  0331H.  A.  Chaikin 
NAVSSC 

1 SEC  610SB.  V.  Dietrich 
I SEC  6105B.  P.  Aaklowitz 
1 SEC  6114B2.  A.  Fuller 
1 SEC  6136.  a.  Eeene 
1 sec  6136.  V.  Fuller 
1 SEC  6136.  D.  Billin«aley 
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CENTE3  DISTRIBUTION 

iTHsaoc 
1 Co<ie  15 
1 Code  1504 
1 Code  1506 
1 Code  1507 

1 Code  1513  - Dr.  Manley  St.  Denis 
1 Code  1524 
1 Code  1552 
1 Code  156 
1 Code  1561 
1 Code  1568 
1 Code  1572 

1 Code  18  > G.  Gleisaner 
1 Code  185  - T.  Corin 
1 Code  1851  * J.  Brainin 
4 Code  1853  - a.  Crons 
1 Code  1853  - B.  ThoMson 
1 Code  184  - U.  Lu8t 

1 Code  1809.3  - D.  Harris 
1 Codie  522 


OTNMOC  IMUCS  TMMEi  TYKS  Cf  NiPOATS 

(1(  OTM8NOC  MVOftTIk.  A POMMAt  SEMltS  PUtUSMiaiG  tNPOMlATlOM  OP 
PUUANEMT  TECNmCAL  VAUJ&.  OiSlGMATSO  IV  A lEIIIAI.  NCfOIIT  NUMKR 

(2t  oepahtmemtal  mwoats.  a SEAUPONMAL  series,  recoroimc  MIPOAMA 
TION  OP  A PRCUMINANV  OM  TEMPOAANV  MATUAC.  OA  OP  LMUTEO  MTEAEST  OA 
SKUMPtCAMCC.  CAAAVIMG  A OVAATMEATAL  ALPHANUMSAC  iOENTiPiCATiOfi 

IS  TECNMCAL  mMCtAAMOA.  AN  MiPOAIIAL  SEAIES.  USUALLY  MTEANAL 
INOAKING  PAPEAt  ON  OfiAfiCT  AEPOATS  TO  SPQNIOAS.  MUHAEAfO  AS  TM  SEAlES 
AEPOHTS.  NOT  POA  GEMgAAL  OISTANHiTiON 


